***********
* Table 1 *
***********

reg agree01 female ageyr edu01  c.ideo01 baseline01 c.scr01 c.inparty01 i.proposal, vce(cluster subject)
reg agree01 female ageyr edu01  c.ideo01 baseline01 c.scr01 c.symp01  i.proposal, vce(cluster subject)
reg agree01 female ageyr edu01  c.ideo01 baseline01 c.scr01##c.inparty01 i.proposal, vce(cluster subject)
reg agree01 female ageyr edu01  c.ideo01 baseline01 c.scr01##c.symp01 i.proposal, vce(cluster subject)

* Robustness check: Ordered Probit Models *
*******************************************

oprobit agree01 female ageyr edu01  c.ideo01 baseline01 c.inparty01 c.scr01 IAparty  i.proposal, vce(cluster subject)
test IAparty

oprobit agree01 female ageyr edu01  c.ideo01 baseline01 c.symp01 c.scr01 IAsym i.proposal, vce(cluster subject)
test IAsym

************
* Fig. 1 *
************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.inparty01##c.scr01 i.proposal, vce(cluster subject)
margins, dydx(inparty01) at(scr01=(0 .2 .4 .6 .8 1)) level (95)
marginsplot, recastci(rline) ciopts(lpattern(shortdash)) plotopts(msymbol(p)) ci1opts(lc(black)) plot1opts(lc(black)) yline(0) ///
title("A. Affiliation as Identification") ///
yscale(range(-.3 .7)) ///
ylabel(-.3(.1).7,nogrid) ///     	
ytitle("Marginal Effect of Identification on Agreement with Party") ///
xtitle("Skin Conductance Response towards Party", height(6)) /// 
	scheme(s2mono) graphregion(lcolor(black) fcolor(white))  	 
graph display, ysize(3) xsize(2.5)

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.symp01##c.scr01 i.proposal, vce(cluster subject)
margins, dydx(symp01) at(scr01=(0 .2 .4 .6 .8 1)) level (95)
marginsplot, recastci(rline) ciopts(lpattern(shortdash)) plotopts(msymbol(p)) ci1opts(lc(black)) plot1opts(lc(black)) yline(0) ///
title("B. Affiliation as Sympathy") ///
yscale(range(-.3 .7)) ///
ylabel(-.3(.1).7,nogrid) ///     	
ytitle("Marginal Effect of Sympathy on Agreement with Party") ///
xtitle("Skin Conductance Response towards Party", height(6)) /// 
	scheme(s2mono) graphregion(lcolor(black) fcolor(white))  	 
graph display, ysize(3) xsize(2.5)

************
* S1 Table *
************

reg agree01 female ageyr edu01 c.ideo01 baseline01 scrv01 scrs01 partyids partyidv i.proposal if imagelist==3, vce(cluster subject)
reg agree01 female ageyr edu01 c.ideo01 baseline01 scrv01 scrs01 symps01 sympv01 i.proposal if imagelist==3, vce(cluster subject)
reg agree01 female ageyr edu01 c.ideo01 baseline01 c.scrv01##c.partyidv c.scrs01##c.partyids i.proposal if imagelist==3, vce(cluster subject)
reg agree01 female ageyr edu01 c.ideo01 baseline01 c.scrv01##c.sympv01 c.scrs01##c.symps01 i.proposal if imagelist==3, vce(cluster subject)

************
* S4 Table *
************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.symp01 c.scr01 IAsym i.proposal, vce(cluster subject)

* Linearity *
*************

acprplot IAsym, lowess

* Outliers *
************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.symp01 c.scr01 IAsym i.proposal
predict tempcooksd, cooksd
display 4/e(N)
summarize tempcooksd if tempcooksd>(4/e(N))

* Robustness check: Exclusion of outliers *
*******************************************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.symp01 c.scr01 IAsym i.proposal if tempcooksd<(4/e(N)), vce(cluster subject)
test IAsym

* Normally distributed error terms *
************************************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.symp01 c.scr01 IAsym i.proposal, vce(cluster subject)
predict tempresid, residuals
qnorm tempresid

* Homoscedascity *
******************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.symp01 c.scr01 IAsym i.proposal
estat imtest, white

* Multicollinearity *
*********************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.symp01 c.scr01 IAsym i.proposal, vce(cluster subject)
estat vif


************
* S5 Table *
************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.inparty c.scr01 IAparty i.proposal, vce(cluster subject)

* Linearity *
*************

acprplot IAparty, lowess

* Outliers *
************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.inparty c.scr01 IAparty i.proposal
predict tempcooksd2, cooksd
display 4/e(N)
summarize tempcooksd2 if tempcooksd2>(4/e(N))

* Robustness check: Exclusion of outliers *
*******************************************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.inparty c.scr01 IAparty i.proposal if tempcooksd2<(4/e(N)),  vce(cluster subject)
test IAparty

* Normally distributed error terms *
************************************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.inparty c.scr01 IAparty i.proposal, vce(cluster subject)
predict tempresid2, residuals
qnorm tempresid2

* Homoscedascity *
******************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.inparty c.scr01 IAparty i.proposal
estat imtest, white

* Multicollinearity *
*********************

reg agree01 female ageyr edu01 c.ideo01 baseline01 c.inparty c.scr01 IAparty i.proposal, vce(cluster subject)
estat vif

*************
* S2 Fig. *
*************

ta partyid

*************
* S3 Fig. *
*************

ta symps01
ta sympv01



